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Abstract 

Among the multiple 5D thick braneworld models that have been proposed in 
the last years in order to address several open problems in modern physics, there is 
one that involves a tachyonic bulk scalar field. In this framework we study a thick 
braneworld with a cosmological background induced on the brane. It turns out that 
the field equations derived from the model with a warped 5D geometry are highly 
non-linear equations that admit a non-trivial solution for the warp factor and the 
tachyon scalar field with a de Sitter 4D cosmological background. Moreover, the 
non-linear tachyonic scalar field, that generates the brane in complicity with warped 
gravity, has the form of a kink-like configuration. However, the non-linear field 
equations' restricting character does not allow one to easily find thick brane solutions 
with a decaying warp factor which leads to the localization of 4D gravity and other 
matter fields. We manage to obtain such a thick brane configuration in this tachyon- 
gravity setup. When analyzing the spectrum of gravity fluctuations in the transverse 
traceless sector we show that 4D gravity is localized due to the presence of a single 
zero mode bound state separated by a continuum of massive Kaluza-Klein (KK) 
modes by a mass gap. This is in contrast with previous results where there is a KK 
massive bound excitation that has no clear physical interpretation. The mass gap is 
determined by the scale of the metric parameter b. We finally compute the corrections 
to Newton's law in this model and show that they exponentially decay and coincide 
with corrections reported in previous results within similar braneworlds with induced 
4D de Sitter metric. 
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1 Introduction 



Within the framework of the braneworld models embedded in a spacetime with extra di- 
mensions and after the success of the thin brane models, where singularities are present 
at the position of the branes, in solving the mass hierarchy and 4D gravity localization 
problems [1, 2], it has became a matter of interest to find smooth braneworld solutions 
(for an interesting review see [3] and references therein). In some models, such solutions 
are obtained by introducing one or several scalar fields in the bulk and the large variety 
of scalar fields that can be used to generate these models gives rise to different scenarios 
[4]-[12]. By following this direction and by using the freedom one has to choose a scalar 
field, several authors have chosen a tachyonic scalar field in the bulk [11]— [15] and addressed 
issues like the mass hierarchy problem, and localization of gravity and matter fields in such 
a model with thin and thick branes. In the original Randall-Sundrum model, a Standard 
Model or TeV brane is introduced at a certain fixed distance, say r c , from the gravitational 
or Planck brane in order to achieve the desired warping, and hence, resolve the hierarchy 
problem in a completely 5D geometrical way. However, this resolution mechanism gives rise 
to a new fine-tuning on the probe brane position, and therefore, to the need of stabilizing 
this brane separation. The stabilization of this brane separation is achieved through the 
Goldberger-Wise mechanism by associating to it a radion scalar field that models the radius 
of the fifth dimension when one ignores the back-reaction of the brane itself [16] . However, 
in [13] it was proved that when one takes into account the full back-reacted system, the 
brane separation becomes unstable when modeled by a standard scalar field even with an 
arbitrary self-interaction potential. 

In order to solve this problem these authors proposed a tachyonic scalar field action for 
modeling and stabilizing the brane separation. They also obtained the desired warping from 
the Planck scale to the TeV one, thus resolving the fine tuning problem of the Higgs mass 
in a stable braneworld scenario and generalizing in a relevant way the Randall-Sundrum 
model (see also [15]). This fact physically motivates the use of a tachyonic scalar field 
within the braneworld paradigm since the back-reaction of the radion field must be taken 
into account in a self-consistent system. 

On the other side, the thick brane configuration constructed in [14] possesses an increas- 
ing warp factor, since most of the attempts to solve the highly non-linear field equations 
leads to imaginary tachyon field configurations. This fact translates into derealization of 
gravity and other kinds of matter like scalar and vector fields, while giving rise to local- 
ization of fermions. Alternatively, the localization of fermionic fields on these thick branes 
was investigated in [17] by using the thick brane solution found in [14] with an auxiliary 
scalar field that couples to the fermionic mass term. In this paper we propose a new thick 
brane tachyonic solution with a decaying warp factor that enables localization of 4D gravity 
as well as other matter fields. By analyzing the dynamics of the metric perturbations we 
see that their spectrum contains a single bound state corresponding to the 4D massless 
graviton of the model and a continuum of KK excitations separated by a mass gap. We 
further compute the corresponding correction to Newton's law by analyzing the influence 
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of these massive modes on the gravitational potential acting between two point massive 
particles located along the center of the thick brane. 



2 The thick brane model and its solution 

The complete action for the tachyonic braneworld model is expressed as follows 
S = J d 5 x^ (J^R - A 5 ) + J d 5 x^V(T)^l + g AB d A Td B T 



(1) 



where the first term describes 5D gravity with a bulk cosmological constant A 5 , the second 
is the action of the matter in the bulk, k 5 is the five-dimensional gravitational coupling 
constant, and A, B = 0,1,2,3,5. Thus, the tachyon field T represents the matter in the 
5D bulk and depends, for simplicity, only on the extra dimension, while V(T) denotes its 
self-interaction potential [18]. 

The Einstein equations with a cosmological constant in five dimensions are given by 

G AB = -Kqab + 4 Tab- (2) 

For the background metric we use the ansatz of a warped 5D line element with an induced 
3-brane with spatially flat cosmological background that reads 



-dt 2 + a 2 it) (dx 2 + dy 2 + cb 2 )] + da 



(3) 



where f(a) is the warp factor and a(t) is the scale factor of the brane. 

The matter field equation is obtained by variation of the action with respect to the 
tachyon and is expressed in the following form 



T + 4/ T (1 + T ) = (1 + T 2 ) 



, 2 ,d T V(T) 



V(T) 



(4) 



By using the ansatz (3) we can obtain the components of the Einstein tensor 

Goo = 3^-3e^(/" + 2/ 2 ), 

Gaa = -2aa-d 2 + 3a 2 e^(/" + 2/' 2 ), 



G n 



-3e 



-2/ 



a 2 



- + ^]+6/ 2 , 



a cr 



(5) 



where "/" and "• " are the derivative with respect to the extra dimension and time, re- 
spectively, while a labels the spatial dimensions x, y and z. The components of the stress 
energy tensor do not depend explicitly on time and read 



T 



bulk 
AB 



g AB V(T) Jl + (VT) 2 - V{T) d A T d B T 

1 + (VT) 2 



(6) 
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while the Einstein equations (2) can be rewritten in a simple way: 



f - r? V{T)T ' 2 c-y ° m 

Thus, we end up with the case of a single brane configuration in a 4D spatially flat cosmo- 
logical background given by equations (4), (7) and (8). 

Now we take the restriction equation (8) and derive it with respect to the extra dimension 
to get 

„ = 4 T'd T V(T) _ 4 T'T"V{T) _ ( a a*\ 

12/V1 + T' 2 12f(l + T' 2 )l 2 [a a 2 J U 
and then equate relations (7) and (9) in order to obtain the following expression 

T" + 4f'T'(l + T' 2 ) + 6//(1+T,2)§ (gJLo) e -2/ = (1 + T ' 2) ^(T) 

which must coincide by mathematical consistency with equation (4). This fact yields an 
expression of de Sitter type for the scale factor since (4) and (10) coincide if and only if 
a(t) = ke bt , where b and k are integration constants. Thus, we must have a de Sitter 4D 
cosmological background defined by 

a(t) = e b \ (11) 

since the constant k can be absorbed into a coordinate redefinition. The last asseveration 
makes clear the role of the action for the tachyonic scalar field as a non-trivial 5D config- 
uration which leads to a braneworld in which the induced metric on the brane is described 
by dSi geometry. 

Here it is convenient to go to conformal coordinates through dw = e~ f ^da which leads 
to the following metric 

ds 2 = e 2f ^ [-dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) + dw 2 ] . (12) 

According to the authors of [14], after obtaining equations (4), (7) and (8), it is easy 
to obtain separated equations for the scalar field T and for the potential V(T) by the 
following procedure. It turns out that despite the high non-linearity of these field equations, 
the derivative of the tachyonic scalar field as well as the arbitrary potential V(T') can 
be expressed in terms of the warp and scale factors of the metric (and their respective 
derivatives) after some simple manipulations: 



r = ±ef 



1 1 (13) 



\ 2 (f 2 + f e 2 f - b 2 
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V(T) 



-2/ 



2 (jn + p + vk e 2/ _ & 2 



/ ,2 + 



A 
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5 „2/ 

e J 



(14) 



where now "/" 



(f 2 + fe 2 f- b 2 ) 

stands for derivatives with respect to w and we have taken into account the 
form of the scale factor corresponding to (11). Thus, by determining a desired behaviour of 
the geometry we completely fix the dynamics of the tachyon field and vice-versa. However, 
the resulting solution must be real and have physical sense. This restriction turns out to be 
very demanding, since several warp factors with "convenient" behaviour lead to a complex 
tachyonic field T and/or self-interaction potential V(T). 
It is straightforward to see that the following warp factor 



/H = ~2 ln 



cosh [b(2w + c)] 



(15) 



where b, c and s > are constants, is a solution of the Einstein and field equations if the 
tachyon scalar field adopts the form 



T(w) = ±, 



2A, 



arctanh 



sinh 



b(2w+c) 
2 



cosh [ b (2w + c)] 



(16) 



while the tachyon potential is given by the following expression 



V(T) 



As 

4 



sech 



-A, T 



\ 



6 sech 



-A, T - 1 



A, 



V24 



1 + sech [b (2w + c)]J2 + 3 sech [b (2w + c)]. 



17) 



In the last two equations we have set s = — 66 2 /A 5 with a negative bulk cosmological 
constant A 5 < 0, a fact that tells us that we have a dS± brane embedded into an AdS 5 [19] 
spacetime. 

Thus, the set of equations (11), (15), (16) and (17) are the necessary and sufficient 
ingredients that give rise to a tachyonic thick brane solution. The warp factor has a decaying 
and vanishing asymptotic behaviour, whereas the tachyon scalar is real and possesses a kink- 
like profile. In contrast with solutions find in [14], here we manage to express the tachyon 
potential V(T) in terms of the tachyonic scalar field T. This potential has a minimum at 
the position of the brane, but it is negative definite as it can be seen from (17). In fact, since 
the tachyon field is bounded, the potential does not extend along the whole fifth dimension 
and remains real as well as bounded. 



3 Gravity localization and corrections to Newton's law 

In order to study the metric and field fluctuations of our system we must perturb both of 
these entities. Since the relevant geometry of the four-dimensional background is de Sitter, 



5 



the fluctuations of the metric may be classified into tensorial, vectorial and scalar modes 
with respect to the transformations associated to the symmetry group dS^ [20]. This fact 
turns out to be very important since at first (linear) order these modes evolve independently 
and, hence, their dynamical equations decouple from each other, even when the perturbed 
Einstein and tachyon field equations were highly coupled non-linear equations. Making use 
of this fact, we shall study the dynamics of the tensorial metric fluctuations which have the 
physical interpretation of the five-dimensional braneworld graviton. These tensorial metric 
fluctuations are gauge invariant and will enable us to determine whether localization of 
four-dimensional gravity, and hence the physics of our four-dimensional world, is feasible 
or not within this model. 

To begin with, we shall write the Einstein equations in the form: 

2 

Rab = ~^59ab + K b T AB . (18) 

where we have made use of the reduced energy-momentum tensor 

V 



T A B = --VJl + iyTfg AB + 



9AB /vt ^2 



(VT) 2 - V A TV B T 



(19) 



In the language of the conformal metric coordinates (12) the 00-component of (18) reads 

/" + 3f> 2 = 3H 2 - 2 -A 5 e^ + fe 2 ^ 2 + (V3T V{T) (2Q) 
3 3 ^1 + (VT) 2 

since 

R 00 = -3H 2 + f" + 3f' 2 (21) 

and 

f V 2f 2 + (VT) 2 

Joo = — e -, =■ (22) 

3 yfl + (VT) 2 

Let us perturb the metric as follows 

ds 2 p = (g AB + h AB )dx A dx B = (g^ + h^dafdx" + e 2f{w) dw\ (23) 

where g AB = (g^e 2 ™) = e 2 ^g AB = e 2 ^)(^,l), = diag[(-l, 5 ija 2 )}, n, v = 
0, . . . , 3 and i, j = 1,2, 3. In this approach we shall impose the transverse traceless condition 
on these metric fluctuations: V^h^ — W — 0, where V is the covariant derivative operator 
with respect to the metric g AB . 

After some algebraic work, the linearized Einstein equations for the transverse traceless 
fluctuation modes adopt the following form 

5Rab = \h AB - ^ + (VT)2 h AB . (24) 
3 3 h + (VT) 2 
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At this stage, it is extremely useful to define a new fluctuation variable = e 2 ^h^ u . By 
making use of the conformal transformation for the Ricci tensor, in terms of the barred 
fluctuations we have 

SR^ = 5R, U - ^j% v - + 3/' 2 ) (25) 

where all the barred quantities are computed with respect to the metric <7ab- On the other 
side, this quantity can also be computed as follows (see, for instance, [21]) 

5R^ = \ (-OV + V A VX + VaV^J - V„V„7^) , (26) 

where Uh^ = g al3 V a V^ u + # 55 V 5 V 5 V- 

In order to compute the involved quantities in this expression we recall that the non-null 
Christoffel symbols are 

T° = a 2 H5 lv rj = H8). 

5 b 

Since T AB = T 5A = and h§A = 0, we have 

V 5 VsV = C> T^cNpKv = n 4 V, (27) 

where n 4 is the d'Alembert operator on dS±. While it is evident that h — implies 
h — from their conformal relationship, it is not easy to see that V^h^ = follows from 
V M /i MI , = and we shall show here that this is indeed the case. On one hand, by definition 
of covariant derivative we have 

= fVaV = 9^{d a v - r^h Ev - r^h Elx ). 

On the other hand, from the relationship between the Christoffel symbols of two conformal 
metrics we get for our case 



^ 7 = and, t 1h-i < ijirii \ /.';;,, 



From both of these relations we conclude that = r*L and, therefore, that V^hn U = 



Since in a curved spacetime the double covariant derivatives do not commute, for the 
the metric fluctuations we have (see [22], for instance) 

VcV^Tf = V,Vch° u + ReX ~ R E ucX (28) 

By making use of the transverse conditions V M /i^ = and the fact that Y AB = Tf A = 0, 
we have 

V c 7? = V 5 hl + VX = 0, V„ VJk = 0. (29) 
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Since h 5 A = and R$abc = R$a = 0, we get 

V c V M /f = R a vK - ifipjit- (30) 
Moreover, for a dS± spacetime the following relations hold 

Rfiva/3 = H 2 (g m g ul3 - g^g va ), R^v = 3H 2 g iiV . (31) 

Hence, we get the following result for the linear variation of the Ricci tensor 

SR^u = g (~n 4 V - V + 8 # 2 V-) ( 32 ) 
Finally, by making use of the relations (20), (24), (25) and the last expression we have 

*4 V + /V + 3f /V - 2 # 2 V = 0, (33) 

under the imposed transverse and traceless conditions V M /i MJ , = 0, h a — 0. 

We further perform the following separation of variables for the metric fluctuations 
h^y = e~2f( w ^( K w)(f) llu (x) that brings the equation (33) into a Schrodinger-like one along 
the extra dimension: 

(~d 2 z + V QM - m 2 ) *(w) = 0, (34) 
where the analog quantum mechanical potential Vqm reads 

V Q M = \f 2 + \f". (35) 

On the other side, the 4D equation that we get from (33) is 

(-d 2 - 3Hd t + e~ 2Ht V 2 - 2H 2 ) <j>{x) = -m 2 0(x), (36) 

where m 2 represents the mass that a 4D observer sees in a de Sitter spacetime [23]- [24] and 
we have omitted the indices of the function ^^{w) for convenience. 

By substituting the expression for the warp factor (15) into (35) the analog quantum 
mechanical potential adopts the form of a modified Poschl-Teller potential that reads 

3b 2 r i 
Vqm = — [ 3 - 7sech 2 (2&u>)J . (37) 

The fact that this potential possesses a definite positive asymptotic value ensures the ex- 
istence of a mass gap in the graviton spectrum of KK massive fluctuations determined by 
|6 2 , or equivalents m = f [24]-[29]. 

By performing the following rescaling of the fifth coordinate v = 2bw we recast (34) into 

- »(«) =(~1) *(»), (38) 



which can be directly compared to the canonical form of the classical eigenvalue problem 
for the Schrodinger equation with a modified Poschl-Teller potential 

-d 2 v - n(n + l)sech 2 w] = E (39) 

with n = 3/4 and £ = fp- — ^- Since ra < 1 there is just one bound state in the mass 
spectrum of KK metric fluctuations: the zero mode massless state which accounts for the 
massless 4D graviton. To the best of our knowledge, this is the first model that presents 
just a single bound state in this kind of spectra within the framework of thick braneworlds. 
Usually one encounters a second bound state which represents a massive KK excitation 
with no clear physical interpretation (see [5], [24]-[29], for instance). The equation (38) can 
be integrated for arbitrary mass and possess the following general solution 

*(«;) = d P 3 M /4 (tanh(26w)) + C 2 Q£ /4 (tanh(2fcw)) (40) 
where C\ and C 2 are constants, while P^ 4 and Q3/4 are associated Legendre functions of 



first and second kind, respectively, with degree v — 3/4 and order p, — yjQ — ^s- For the 
massless case the zero mode has p, — v — 3/4 and the exact solution can be expressed as 



P^ (tanh(26w)) + - Q^J (tanh(26w)) 



(41) 



where now fci = — C\ > and we have set C 2 = 2Ci/n in order to get a localized configura- 
tion. This bound state is physically interpreted as a stable graviton localized on the brane 
since there are no states with negative squared masses due to the structure of the potential 
(39). Finally, there is also a continuum of KK massive modes in the spectrum that starts 
from m > 36/2 and is described by eigenfunctions with imaginary order ±/i = ±ip [26]: 

* m (w) = £C ± P±f (tanh(2H) , (42) 
± 



where C± are arbitrary constants and p = y fp- — J6- These KK massive modes must behave 
as plane waves asymptotically [24, 25]. This fact can be seen by considering masses 2m > 3b 
and taking into account that the constants C± depend on p in general, i.e. 

= C ± {p)P^ (tanh(2H) • (43) 

We further make an expansion for large w of the argument of the associated Legendre 
functions of first kind: 

tanh(26w) ~ 1 - 2e~ 4bw (44) 

and compute the asymptotic behaviour of P^Jl (tanh(26w)) according to formula (8) of 
section 3.9.2 of [30]: 

Ptjl (tanh(2H) ~ (45) 
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The normalization condition of (43) in the plane wave sense leads to the following normal- 
ization constants 



C+(p) = C_(p) 



|r(i + v)| 



2tt 



(46) 



since |T(1 — ip)\ = |T(1 +ip)\. Thus, by substituting (45) and (46) into (43) we obtain the 
asymptotic behaviour of these associated Legendre functions: 



1 



D ±i2bpw 



2tt 



(47) 



which corresponds to plane waves as expected. 

Once we have at hand an analytical expression for the KK massive modes, we should 
be able to compute the corresponding small corrections to Newton's law due to these 5D 
massive modes. This is achieved by taking the thin brane limit b — > oo, locating a 
probe mass, Mi, in the center of the brane in the transverse direction and computing the 
gravitational potential generated by this particle that another massive particle with M 2 
realizes. The corrections to the Newtonian potential generated by massive gravitons in the 
thin brane limit can be expressed as follows [6] 



U(r) 



MiM 2 



G 4 + M, 



L 



oo 

' dme~ 

'mo 



MiM 2 



(G 4 + AG 4 ) , (48) 



where w = w sets the position where the brane is located, m = 36/2 for our case, G 4 is 
the gravitational 4D coupling constant and ^^(wq) denotes the continuum of KK massive 
modes that must be integrated over their masses in order to get the searched corrections. 
Now we proceed to calculate |^(0)| 2 at w = and get 



|^(0)| 



V(l+ip 



r(¥+f)r(|+f) 



(49) 



where we set v = 3/4 and made use of the following relation for the associated Legendre 
functions of first kind [31]: 



P*(0) = 



2^ 



2 



(50) 



r(±^)r(i 

Substitution of (49) into the second term of (48) yields the following expression for AG 4 : 

V(l + ip 



AG 



PCX) 

4 = M~ 3 / dm e 

J mo 



rftKf)r(f+f) 



(51) 
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In order to calculate this integral it is convenient to perform a change of variable from m 
to p: 



AGa = 2M~ 3 b 



f 

Jo 



9 

16p 2 



r(i + i P ) 



r (¥+f) r (I+f) 



dp. 



(52) 



This integral can be calculated in the thin brane limit b — > oo, where it is dominated by 
the region of small - p. Thus (52) can be well-approximated by expanding the factor that 
multiplies the exponential at p = [26]. Thus, the contribution to Newton's law made by 
the continuum of KK massive modes in the thin brane limit reads: 



AGa 



M. 



-3 



e 2 



br 



l + O 



(53) 



These corrections are exponentially suppressed as in other braneworld models with an 
induced 4D Minkowski [32], (see also [26]-[27]) or de Sitter metric [28, 29]. 



4 Discussion 

In this paper we have presented a thick braneworld model generated by a tachyon scalar 
field coupled to gravity with a bulk cosmological constant and a de Sitter metric induced on 
the brane. We were able to obtain an exact solution with a decaying warp factor that gives 
rise to localization of 4D gravity when studying the metric fluctuations. The structure of the 
corresponding graviton spectrum is novel in the sense that it contains just one bound state 
(the massless zero mode) which is physically interpreted as the 4D graviton, separated by 
a mass gap from a continuum of KK massive excitations. We got an explicit expression for 
these KK massive modes, a fact that enables us to analytically compute the corrections to 
Newton's law coming from the extra dimension. As we mentioned above, these corrections 
are exponentially suppressed and coincide with previous results reported in the literature 
within the framework of similar braneworld models with the same 4D de Sitter induced 
metric. 

Since we are considering a braneworld with an induced 4D de Sitter metric, in principle, 
we are able to reproduce the early inflation and accelerated expansion epochs of our universe 
within our model. However, a more general and realistic ansatz for our setup that attempts 
to describe the late time behaviour of our 4D universe should involve both a time-depending 
tachyon field and a time-depending warp factor, since from the cosmological viewpoint one 
needs to obtain scale factors that reproduce in a better way (closer to the observations) the 
accelerated expansion of the universe (which takes into account its dark matter component) 
than the de Sitter metric does. This line of research is under current investigation and could 
lead to interesting new results in cosmology. 
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